Enhanced Cosmological GRB Rates and Implications for Cosmogenic Neutrinos 
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biased tracer of the SFH B221I23I1 . In particular, the host galax- 
ies of GRBs have a distinct tendency to be subluminous lf22ll 
and metal-poor [24]. This has been demonstrated for GRB 
hosts both locally M25ll and at cosmological distances j26ll . 
which suggests that low metallicity is a key ingredient in the 
production of a gamma-ray burst 112511 . As we discuss in Sec- 
tion [II] a rapidly rotating star, as required in the collapsar 
model [27], can retain much of its original mass and angu- 
lar momentum if it is metal-poor |28]. An anti-correlation 
with metallicity would im ply that the cosmological GRB rate 
peaks at a higher redshift [29, 30], which now seems to be in- 
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I. INTRODUCTION 

The origin of ultrahigh energy cosmic rays (UHECR) is 
one of the great remaining mysteries in astrophysics [ lj. The 
cosmic-ray spectrum has been measured to beyond 10 19 eV, 
with a number of events with energy exceeding 1 20 e V 0, Ht] ; 
however, it is still debated how such highly energetic particles 
can be produced. It is now generally accepted that UHECR 
are of an extragalactic origin 0]. However, above ~ 10 195 eV, 
the process of photopion production (p 7 — > Nn) on the cos- 
mic microwave background (CMB) is expected to lead to a 
significant diminution of the cosmic-ray spectrum, the well- 
known GZK effect HI] . The relatively short attenuation length 
associated with the GZK process (3 0] necessitates that the 
observed UHECR arise from local sources. 

The decay of charged pions produced in this process re- 
sults in a flux of ultrahigh energy cosmogenic neutrinos JHIgl]. 
While the observed UHECR spectrum is somewhat insensi- 
tive to variations in cosmic source evolution 1 10], the cosmo- 
genic neutrino flux can be greatly enhanced by strong evo- 
lution with redshift |9l Ell, as neutrinos can be produced 
in larger quantities due to the decreased photopion threshold 
(since 7cmb x 1 +Z), and can themselves traverse cosmologi- 
cal distances without attenuation. 

Few classes of astrophysical objects can possibly account 
for the observed cosmic-ray spectrum 0|. Active galactic nu- 
clei (AGN) have long been considered as possible UHECR 
sources 1131 . Relatively recently, a potential connection be- 
tween gam ma-ray bursts (GRBs) and UHECR has been ex- 
plored H13L fl4l fl5ll . GRBs, which are now generally accepted 
to be related to core-collapse supernovae flql . are violent 
events which release great amounts of energy in the form of 
gamma-rays (~ 10 51 — 10 52 erg) |[T7ll . A number of models 
have been proposed to utilize their ultra-relativistic environ- 
ment to accelerate protons to energies > 10 20 eV|E[ll[l!. 
It has also been noted that the source emissivity required to 
account for > 10 19 eV cosmic rays is comparable to that of 
gamma-ray bursts lfl3l,[l8ll . 

Just as the core-collapse supernova rate density seems to 
follow the cosmic star formation history (SFH) Il20ll . the same 
might be expected of the cosmological GRB rate density ll2"lll . 
However, there is mounting evidence that GRBs are not an un- 
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FIG. 1: Models of cosmic-ray source evolution (i.e. yield vs. z, nor- 
malized to 1 at z = 0). From top-to-bottom, the metallicity-dependent 
GRB rate density (this work), the quasar (QSO) evolution model used 
in Refs. iTll u%, and the fit to the cosmic star formation history 
(SFH) of Ref. 1 20]. Models similar to the latter two have been fre- 
quently used in cosmic-ray studies. 
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dicated by Swift observations (3lll32Tl . Simply put: the metal- 
licity of the universe decreases with redshift, which implies a 
stronger evolution of the GRB rate density than would be ex- 
pected from the SFH alone. 

When GRBs are considered as the source of UHECR (with 
identical cosmic-ray production per burst), the change in the 
cosmological cosmic-ray emissivity is simply determined by 
the burst rate history. In the metallicity-biased GRB model, 
this evolution is quite strong, as illustrated in Fig. Q] even 
exceeding that of models used in past cosmic-ray studies, 
which have traced quasar (QSO) luminosities H33I1 or the 
SFH. We examine the effect of enhanced GRB rate evolu- 
tion on the expected flux of cosmogenic neutrinos, the mea- 
surement of which may provide the only way to break de- 
generacies between cosmic -ray models [34]. In addition to 
naturally explaining the abundance of high-redshift bursts ob- 
served by Swift, this strong evolution leads to a measurable 
neutrino signal, improving the near-term prospects for assess- 
ing this scenario with upcoming ultrahigh energy neutrino 
detectors ll35l l36l l37tl . This result would still hold, in gen- 
eral, even if another mechanism is ultimately shown to ac- 
count for the Swift results through enhanced GRB evolution. 
Though we focus on the impact on cosmogenic neutrino pro- 
duction (as Dermer and Holmes have also recently done for 
GRB rates directly tracing several SFHs [38]), the metallicity- 
biased evolution model would also have implications for pre- 
dictions of diffuse neutrino fluxes directly produced in GRBs 
(e.g., Refs. idllll]) and their prospects for detection j40ll4lll . 



II. METALS AND THE PREDICTED GRB RATE 

Any treatment of cosmological cosmic-ray or neutrino pro- 
duction must account for source evolution. Traditionally, 
the most commonly used cosmic-ray evolution models have 
tracked quasars [ 33] or the SFH 112011 (particularly in GRB- 
related studies). Observations indicate, however, that GRBs 
do not faithfully trace the SFH [22, 23]. In fact, it now appears 
that the GRB rate density is actually evolving more strongly 
than the standard SFH ll3ll[32ll42Tl . A natural explanation may 
be found by considering the properties of GRB host gala xies 
in the context of the single-star collapsar model j27l l28ll . In 
this model, GRB progenitors are rapidly-rotating Wolf-Rayet 
stars which undergo a core-collapse event that produces a 
black hole (or possibly a rapidly -rotating neutron star) j43ll . 
After collapse, this rapid rotation allows for the formation of 
highly relativistic jets which, when viewed on-axis, are ob- 
served as a burst 114411 . 

Observationally, all known supernova counterparts of 
GRBs are Type Ic, with the implication that the dying star 
lacked an outer hydrogen/helium envelope Olal . The winds 
of Wolf-Rayet stars, which are typically the cause of the loss 
of this envelope, are known to increase in strength with stel- 
lar metallicity (particularly iron) 04511 . Importantly, in this 
wind-induced loss process, angular momentum (which is par- 
ticularly important for forming jets Mm ) is lost along with 
mass 114711 - This loss of angular momentum can be avoided 
if the progenitor has a very low metal content. In addition 
to having weaker winds, a rapidly-rotating, metal-poor mas- 
sive star can avoid the production of an envelope altogether by 
completely mixing its interior, which results in the hydrogen 



being circulated into the core and burned |28|. This would be 
impossible if the star was not metal-poor, as stellar mixing is 
expected to be inhibited by a high metal content j48ll . 

Studies of GRB hosts demonstrate that these galaxies tend 
to have a low metallicity 1251 12611 . Compared to consider- 
ations of galactic luminosity observations alone [49], direct 
spectroscopic measurements reveal metallicities that are sig- 
nificantly lower than expected [25]. While the metallicity of 
the GRB progenitor may not be directly measured, the metal- 
licity of the galaxy itself should be indicative. A connec- 
tion between GRBs and metallicity implies that the cosmo- 
logical GRB rate should be dependent upon the formation 
rate of metal-poor stars. At higher redshifts, the universe is 
less metal-enriched than at the present [50], resulting in an 
increase in the expected evolution of the GRB rate density 
(compared to the SFH) as shown by Langer and Norman 112911 
(see also Refs. lHjHH)]). Accounting for stellar evolution ef- 
fects, Yoon, Langer, and Norman calculated the expected ra- 
tio of GRBs to core-collapse supernovae as a function of red- 
shift, as seen in Fig. 7 of Ref. [30], which we approximate as 
"grb(z) x (1 + z) 1 ' 4 «sn(z). This is consistent with estimates 
of this ratio [53] and can be used to calculate the increase in 
the absolute GRB rate density (of importance to cosmic -ray 
studies). 

As the progenitors of core-collapse supernovae are very 
short lived, «sn(z) is expected to closely follow the SFH (as- 
suming an unchanging IMF O). Utilizing constraints from 
the diffuse supernova neutrino background 15511 (in addition to 
direct observations), an updated SFH was derived by Hopkins 
and Beacom M20I1 . which is well-described by the piecewise- 
linear fit: «sfh(z) x (1 +z) a , where a — (3.4, 0, -7) when 
(z < 1, 1 < Z < 4.5, 4.5 < z). Combining this SFH with the 
parametrized form of the GRB/SN ratio, we find the source 
evolution term, #grb(z), to have the form 

(1+z) 4 - 8 : z< 1 
^grb(z)« (1+z) 14 : Kz<4.5 (1) 
(1+z)- 5 - 6 : 4.5 <z, 

with «grb (z) = «grb (0) x Wgrb (z) • In Fig. [T] we present this 
evolution model, in comparison to the SFH alone, as well as 
the quasar (QSO) evolution model used in Refs. ifllL 11911 . Note 
that our rate history evolves as Wgrb{z) ~ (1 + z) 4 ' 8 for z < 
1, which is substantially steeper than the models commonly 
considered. In particular, the QSO model, which has been 
extensively used, only rises as Wqso{z) ~ (1 +z) 3 - Previous 
studies have found an evolution in GRB luminosity °< (1 + 
z) 114211 . which would result in a similar effect on cosmic- 
ray evolution. While this may just be due to stronger beaming, 
we will not consider this scenario and instead focus only upon 
changes in the GRB rate density. GRB rates that evolve more 
strongly than the standard SFH have, in fact, recently been 
determined using Swift data | [3"Tl l32ll . 

To directly check for compatibility with observations, 
we calculate the expected redshift distribution (assuming 
Swift 1 56] detection sensitivity) for the metallicity-dependent 
evolution, following the GRB model presented by Le and Der- 
mer ll3lll . and compare to a sample of 46 long bursts discov- 
ered by Swift (compiled from the updated lists of Ref. Iplll ). 
all of which have reliable redshifts. Details of this calcula- 
tion are contained in AppendixlAl In particular, we assume an 



3 




FIG. 2: Predicted redshift distributions of observable gamma-ray 
bursts from the strongly-evolving (metallicity-biased) model (solid 
line) and the SFH alone (dotted line), compared to Swift bursts with 
known redshifts. In the top panel, the cumulative burst distribution is 
shown (shaded region). The bottom panel shows the differential dis- 
tributions (with arbitrary normalizations), which allow for a direct 
comparison with the data (in bins of width 0.5). The excess of low-z 
events may be due to a LLGRB class (as discussed in SectionlTIIl. 



these events are typically weaker than cosmological bursts, 
they may be several orders of magnitude more abundant in the 
local universe jrJTIl . One possibility is that these LLGRBs are 
the result of progenitors that had a higher metal content than 
typical GRBs, leading to a less luminous (and less beamed) 
gamma-ray output. In the context of the collapsar model, this 
is not an unreasonable conclusion, as higher progenitor metal- 
licity should lead to increased loss rates of both mass and an- 
gular momentum. 

One consequence of the strong evolution that we present 
is a predicted burst rate density of «grb(z = 0) ~ 4 Gpc~ 3 
yr _1 in the local universe (when normalized to the observed 
redshift distribution), which quickly rises to ~ 20 Gpc~ 3 yr _1 
by redshift z ~ 0.4. The average emissivity from these bursts 
(at z < 0.4) is <r G RB - e r x 10 Gpc~ 3 yr" 1 - 5 x 10 43 erg 
Mpc -3 yr _I , which is comparable to the emissivity required 
to account for the > 10 19 eV UHECR flux, <? C r ~ few x 10 44 
erg Mpc~ 3 yr~' (as found in Ref. fl3ln . Considering that 
there is no a priori reason for these numbers to be so simi- 
lar, along with the fact that the region of the cosmic-ray spec- 
trum of greatest interest (> 10 19 eV) must arise from sources 
at z < 0.4 (as we will discuss in SectionHITl), this result is quite 
interesting. GRBs may also be more luminous in cosmic rays 
due to a baryon loading factor (/cr) that may be > 10 ll62l 16311 . 
which could account for any difference. These considerations, 
along with the isotropy of the measured cosmic-ray spectrum, 
allow for GRBs to be further examined as a candidate to pro- 
duce the observed UHECR. 



opening angle that varies between m ; n = 0.05 and max = 0.5 
(with a power law distribution of the form « (1 — cos0) s 
where the index is s ~ —1.5) and an absolute gamma-ray en- 
ergy per GRB of £y = 5 x 10 51 erg released in ~ 8t = 10 s, all 
of which are well within the reasonable ranges of these param- 
eters 115711 assuming a uniform jet (see Ref. 03811 for alternative 
jet models). 

In Fig. 12 we compare the distribution of Swift bursts to the 
predictions of the metallicity-biased model and the standard 
SFH. The top panel suggests that the cumulative distribution 
of observed bursts (shaded region) can be approximately de- 
scribed by either the metallicity-dependent evolution or the 
SFH. However, another perspective is presented in the bot- 
tom panel, where the differential event distribution is shown 
(in redshift bins of width 0.5). This allows for a more di- 
rect comparison with redshift data, which reveals that most 
bursts are discovered between z ~ 1 — 4, in good agreement 
with our strongly evolving rate. This stronger evolution al- 
lows for a better fit to the data, when compared to the SFH. 
A more detailed discussion of the uncertainties involved in 
GRB detection (particularly the observed bursts that lack a re- 
liable redshift), as well as the various degeneracies that exist 
between parameters in GRB modeling, is beyond the scope 
of this work. We refer the interested reader to Ref. 113111 (and 
references therein). 

While the metallicity-biased model is in generally good 
agreement with Swift observations at high redshifts, an excess 
of events seems to be present at low redshifts, which affects 
the cumulative distribution. This might be explained by a sep- 
arate class of low-luminosity GRBs (LLGRB) [59, 60]. While 



III. UHECR AND GRBS 

To be considered as a viable source of UHECR, gamma- 
ray bursts must have the ability to both accelerate protons to 
energies > 10 20 eV and generate a cosmic-ray flux adequate to 
explain the observed spectrum. In conventional GRB models, 
a portion of the kinetic energy of a relativistically expanding 
fireball (with Lorentz factor F ~ few hundred) is converted 
into internal energy Il64ll . Electrons are accelerated inside this 
jet by internal shocks and subsequently produce gamma-rays 
through synchrotron and inverse-Compton processes ll65ll . 

Protons should also be shock-accelerated in a similar fash- 
ion. In the internal shock model, the shocks that accelerate 
protons are expected to be only mildly relativistic in the wind 
rest frame, resulting in an ~ E~ 2 spectrum. In order to effi- 
ciently accelerate protons to ultrahigh energies, the time scale 
of acceleration should be shorter than both the wind expan- 
sion time (to allow for an adequate period of confinement in 
shocked regions) and the proton energy loss time scale. The 
former sets the ratio of magnetic field and electron energy den- 
sities to order unity, which is necessary in order to account for 
gamma-ray emission from synchrotron emission boosted to 
the observer's frame. The latter imposes an upper limit on 
magnetic field strength (and lower limit on Lorentz factor). 
See Refs. 11131 fl4l [1511 for details (and alternative models). 

Cosmic rays mainly lose energy through three processes 
during propagation. At very high energies (> 3xl0 19 eV), 
UHECR energy loss is dominated by photopion production on 
the CMB, p y — > N %, which has a maximum cross section at 
the A(1232) resonance ll6l l66ll . Below the photopion threshold 
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(and with E p > 10 18 eV), electron-positron pair production 
on diffuse photons, pj —> pe + e~, becomes important l67ll . 
While the cross section for this process is high, each interac- 
tion results in only a small energy loss. Finally, the expansion 
of the universe results in adiabatic energy loss, which is inde- 
pendent of energy. 

We can account for the total proton energy loss through 
the characteristic timescales associated with each process (as 
in Refs. M H): tf\E p ,z) = T K l {E p ,z) + x^(E p ,z) + 
X~ (z). The rate of energy loss for propagating protons is then 
dlnEp/dt = Tj l (E p ,z). With the redshift-dependent energy 
loss rate determined, the injection energies of cosmic rays, 
E' p = E ! p (Ep,z), can be calculated as a function of detected 
energy (E p ) and originating redshift (z) through the differen- 
tial equation 



1 



1 



1 dE„ 



E p dz %t{E p ,z) dz/dt 



(2) 



where dz/dt = H (l +z)p M (l +z) 3 +Q. A } 1/2 (with Q. M = 
0.3, H A = 0.7, and H = 70 (km/s)/Mpc). 

The top panel of Fig. [3] shows cosmic -ray energy loss with 
redshift. The lines correspond to the injection energy needed 
at redshift z in order for a cosmic ray to be observed with 
a given energy at Earth (z = 0). For example, a cosmic-ray 
proton with a measured energy of 10 19 eV must have been 
produced at z < 0.4. This effect is analogous to the Fazio- 
Stecker relation for gamma rays 16911 . As energy losses above 
3xl0 19 eV are quite severe, the observed UHECR must be 
produced in the local universe. This, along with the neces- 
sary conditions required to produce UHECR at all, strongly 
constrains the population of prospective sources. 

In evaluation of the cosmic -ray spectrum, we have assumed 
that, at these energies, the cosmic-ray spectrum is entirely 
composed of protons with an injection spectrum (p(E' p ) = 

jYEp ^ (per unit comoving volume per unit energy, per unit 
time), which is cut-off at a chosen E cut and normalized to ac- 
count for the observed spectrum of UHECR. Cosmic rays in- 
jected at E' t will experience energy losses and be detected at 
E p . Taking into account the evolution of sources, W(z), we 
can calculate the UHECR spectrum as, 



dN p 
dE„ 



c 

An 



<P{E') 



dE p W{z) 
dE p dz/dt 



dz. 



(3) 



We calculate the partial derivative dE' p (E,z)/dE p numeri- 
cally from Eq. [2] In the bottom panel of Fig. [3] we present 
the expected cosmic -ray spectrum for the three source evolu- 
tion models presented in Fig. Q] with an injection spectrum 
of the form E~ 2 (plotted as E 3 x Flux to emphasize spectral 
features ll70ll ). We choose the normalization of (p, JV , such 
that E 3 x dNp/ dEp = 2 x 10 20 eV 2 cm" 2 s" 1 sr" 1 at 10 19 eV 
for these spectra (and in our subsequent results), as the un- 
certainties at higher energies await the resolution that will be 
delivered by Auger [3]. This usually necessitates the local 
emissivity of sources between 10 19 — 10 21 eV to be on the or- 
der of # C R ~ 5 x 10 44 erg Mpc" 3 yr" 1 OS [El]. 

Note that the data from the AGASA, HiRes and Yakutsk 
experiments [2] differs by almost a factor of 2 in overall nor- 
malization in our E 3 plot (light shaded region), however, the 




FIG. 3: Top: Cosmic-ray energy loss with redshift. Lines illustrate 
the injection energy required at redshift z in order to be detected at a 
given energy at z = 0. Shown for illustrative purposes are the realms 
of photopion (dark-shaded region) and pair-production (light-shaded 
region) losses. Bottom: Cosmic-ray spectra expected from the GRB 
(solid), QSO (dashed), and SFH (dotted) source evolution models, 
assuming y = 2 and normalization at E = 10 19 eV. All curves are 
well within range of experimental data above 10 19 eV (light-shaded 
region, dark-shaded when normalized to the spectral dip). 



overall spectra do not significantly disagree [71] and model- 
independent information can be extracted from the shape 
alone, especially when various experimental data are normal- 
ized to the spectral dip feature at ~ 10 19 eV l68ll (dark shaded 
region). 

When considering the cosmic -ray spectrum resulting from 
GRB evolution, we have followed the general form of Wax- 
man by assuming that the entire UHECR spectrum above 



10 19 eV can be accounted for by GRBs, with an £ 



injec- 



tion spectrum lfl3l [l8ll . However, for an injection spectrum of 
this form, it is challenging to account for cosmic rays with en- 
ergies < 10 19 eV iH. Strong source evolution with redshift 
offers a certain degree of assistance. It has been assumed pre- 
viously that an extension of the Galactic cosmic-ray spectrum 
may account for any remaining deficiency 117311 . Contributions 
from extragalactic sources with a lower energy reach may also 
be considered (e.g., AGN, cluster shocks Il74ll . or LLGRBs). 
If LLGRBs are able to produce cosmic rays (as suggested in 
Ref. 117511 ). they may be able to contribute to the flux of lower 
energy cosmic rays (since they have a high local rate jrJTIl ) and 
alleviate the need for an extension of the Galactic component. 
A recent Galactic GRB has also been proposed as a potential 
source of lower energy cosmic rays II76I1 . While the Milky 
Way is generally metal-rich (contrary to observations of typ- 
ical GRB hosts), it possesses a gradient in its metallicity dis- 
tribution which may have allowed for a low-luminosity burst 



in the past. Alternatively, one can consider an E injec- 
tion spectrum (as might be expected from relativistic shocks), 
which may offer a better fit at lower energies 11381 17611 . This 



5 



would require a larger fcR, which itself may explain GRBs 
with low radiative efficiencies 17611 . 

It is important to note that the cosmic -ray spectrum it- 
self depends only mildly upon evolution (since the observed 
UHECR with E > 10 19 eV must be produced at z < 0.4) and 
may be explained by a combination of GRBs and a lower- 
energy component. However, the strong evolution implied by 
the metallicity-biased GRB model results in an increased flux 
of cosmogenic neutrinos. This neutrino flux, which must be 
present if GRBs are to account for the observed UHECR, pro- 
vides an independent test of the GRB-UHECR model, which 
we will now examine. 



IV. IMPLICATIONS FOR COSMOGENIC NEUTRINOS 

The flux of cosmogenic neutrinos produced via the GZK 
process is quite sensitive to cosmic-ray source evolution i34Tl . 
This is due to the unique ability of neutrinos to propagate 
through large distances at very high energies without appre- 
ciable energy loss (unlike cosmic rays), combined with a 
lower threshold for photopion production at larger z (since 
Tcmb 05 1 +z). In fact, a significant portion of the expected 
neutrino flux originates between from z ~ 1 — 4. 

Near the photopion production threshold, about 20% of the 
original proton energy is lost in each interaction. Neutrinos 
are subsequently produced through the decay chain, 7T+ — > 
jl + V M — > e + v e Vji Vfi, with each daughter neutrino (v e ,V^,V^) 
receiving ^1/20 of the parent proton energy [ 77] . At higher 
energies, multi-pion production dominates 16611 : however, this 
approximation remains viable, as the inelasticity of each in- 
teraction increases (approaching 50%). 

Waxman and Bahcall (WB) have presented an upper bound 
on cosmogenic neutrino production (shown in Figs. [4] and [5] 
as a shaded band) based on the assumption of an E~ 2 injec- 
tion spectrum, with normalization chosen between 10 19 eV to 
10 21 eV to produce the observed cosmic-ray spectrum lfl9tl . 
This yields an energy-dependent rate of cosmic -ray genera- 
tion of jKmb ~ 10 44 erg Mpc~ 3 yr _1 (with <f W B ~5x */*wb)- 



yr 

The total + energy flux at Earth (not corrected for oscil- 
lations) is estimated as 



' v dE v 



C A, 1 E 



x 15- 



eV 



(4) 



wheref// rj 10 10 yristhe Hubble time and the factor l/4arises 
from the assumption that only one quarter of the energy lost 
is carried away by muon neutrinos. Adiabatic redshift losses 
and the effect of source evolution are taken into account by % 
(estimated to be | ~ 0.6 with no source evolution, ~ 3 with 
QSO-like evolution) OH]. 

The neutrino spectrum produced in the GZK process may 
be better approximated for a general cosmic -ray injection 
spectrum and source evolution through a somewhat more so- 
phisticated approach. As the energy loss distance at these en- 
ergies is relatively short, we assume that cosmic rays lose all 
of their energy rapidly. The fraction of the original proton en- 
ergy that is lost to neutrinos can then be parametrized with 
a gradual step function, S"{E) = 0.45/(1 + (E t /E) 2 ), where 
0.45 is the asymptotic fraction of injected cosmic -ray energy 
transferred to neutrinos above > 10 21 eV (as shown in Fig. 1 



of Ref. |0) and E, ~ 2 x 10 20 eV governs the onset of photo- 
pion production above ~ 3 x 10 19 eV. The total neutrino flux 
at Earth can be cast as 



dN v 



dE v dz/dt 



dz, 



where E' p 



(5) 

20(1 +z)E v and the factor 20 reflects the ap- 
proximation of each daughter neutrino receiving about 1 /20 
of the injected proton energy. A more detailed derivation of 
this formula is given in Appendix iBl The additional factor of 
(1 + z) in 5? accounts for the lowering of the photopion energy 
threshold as the CMB temperature increases at higher redshift. 
This simple formulation is similar to the notation based on 
neutrino yield functions often used in prior studies |[TTll78lF79ll 
and provides a reasonably accurate neutrino spectrum which 
agrees rather well with the literature (in the energy range most 
interesting to UHE neutrino detectors), with a deviation from 
the simulated spectra not larger than the variations typically 
seen between such simulations. The overall normalization 
(set by <p) is again chosen such that the corresponding pre- 
dicted cosmic -ray flux agrees with the measured cosmic -ray 
data at 10 19 eV (where it is well-determined), as discussed in 
Section Hill 

Given an UHECR source evolution, W(z), we can calculate 
the expected flux of cosmogenic neutrinos produced through 
photopion production on the CMB. While our simple, ana- 
lytic method allows for a more transparent look at the effects 
of source evolution, it does not fully encompass the particle 
physics involved, particularly the low and high energy tails of 
the distributions of particle decays (which affect the low- and 
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FIG. 4: Expected (all-flavor) cosmogenic neutrino fluxes assuming 
various evolution scenarios. From top-to-bottom, are the fluxes re- 
sulting from the strongly evolving (metallicity-dependent) GRB rate 
density, QSO-like evolution, and the SFH. Shown for comparison 
are the Waxman-Bahcall bound (shaded band) and the expected sen- 
sitivities for ANITA (diamonds) and ARIANNA/SalSA (triangles). 



6 



7=1.5,2,2.5 : 
E =10 21 eV 




io- 20 - 

10 10 



10" 



E [eV] 



FIG. 5: Expected (all-flavor) cosmogenic neutrino fluxes result- 
ing from various UHECR injection spectra and assuming strong 
(metallicity-dependent) GRB evolution. From top-to-bottom, are the 
fluxes with spectral indices 1.5, 2, and 2.5. Shown for comparison 
are the Waxman-Bahcall bound (shaded band) and the expected sen- 
sitivities for ANITA (diamonds) and ARIANNA/SalSA (triangles). 



in a manner that would not be possible by cosmic-ray obser- 
vations alone. 

We can also consider the effect of varying the assumed pa- 
rameters in the injection spectrum. Since only cosmic rays 
with energies greater than 10 19 eV will ever contribute sig- 
nificantly (even at high redshifts) to the neutrino flux, harder 
cosmic-ray spectra will result in larger fluxes. Shown in Fig. [5] 
are the neutrino fluxes arising from cosmic-ray injection spec- 
tra of 7 = 1.5, 2.0, and 2.5, assuming the strong evolution. 
In particular, a 7= 2.2 spectrum might arise from relativis- 
tic shocks (see Ref. [38] for details). Since y < 2 would be 
quite difficult to reconcile with cosmic -ray observations, this 
range can be regarded as an upper bound, which will soon 
be constrained by ANITA. Assuming a lower cutoff energy 
would be similar in effect to a softening of the spectral in- 
dex. Another interesting scenario is producing a softer spec- 
tral index (7 > 2), even if each GRB possesses an intrinsic 
E~ 2 spectrum. As proposed in Ref. [83], if the distribution of 
the high-energy cutoffs of cosmic -ray sources follows a power 
law, then the overall spectrum that is observed will follow 
a power law with a different slope. Such an argument may 
be particularly applicable to GRBs, as cosmic rays might be 
produced by a spectrum of sources ranging from LLGRBs to 
strong GRBs with very metal-poor progenitors. 



V. DISCUSSION AND CONCLUSIONS 



high-energy ends of the neutrino spectrum). We utilize the 
publicly -available simulation package CRPropa [80], which 
uses the SOPHIA B81I1 code to handle particle processes, for 
this purpose. We have made use of the analytical estimate de- 
scribed above and an extensive comparison to previous results 
presented in the literature (with similar parameters), to verify 
the results. At the redshifts of greatest interest, z ~ — 4, only 
cosmic rays with £>5x 10 / (1 +z) eV can ever have the 
ability to produce neutrinos through the GZK process (even 
with a decreased photopion threshold). At lower energies, 
photopion production can be facilitated by the cosmic infrared 
background (IRB), resulting in additional neutrinos of cor- 
respondingly lower energy 1831 . Additionally, extragalactic 
magnetic fields may increase the path length of cosmic -ray 
propagation [7]. However, due to the various uncertainties in- 
volved, we have chosen not include these effects. 

Fig. 2] compares the resulting neutrino flux from the GRB 
model to those expected from QSO-like evolution and the 
SFH, assuming an E~ 2 injection spectrum (with a sharp cutoff 
at 10 21 eV). The flux shown is the sum of all neutrino flavors 
(both electron and muon types are produced at the source), as 
a detector such as ANITA has a nearly equal taste for each 
flavor. As can be clearly seen, the stronger evolution of the 
GRB model produces a larger neutrino flux than other mod- 
els. Also shown are the expected sensitivities of ANITA 113511 
(which began taking data in late 2006) and the proposed AR- 
IANNA/SalSA detectors ll36Tl . The flux predicted by strong 
GRB evolution, with parameters given by our metallicity- 
biased model, extends well into the reach of ANITA. This 
presents an opportunity (within the next few years) to either 
confirm or significantly constrain the parameters of this model 



The increase in our understanding of gamma-ray bursts can 
be traced to the improved capabilities now available to study 
these phenomena. The ability to quickly and accurately local- 
ize a GRB has led to the establishment of a GRB-supernova 
connection and allowed for the study of the host galaxies in 
which these events occur. Observations indicate that these 
hosts tend to be underluminous, star-forming and metal-poor. 
A connection between a GRB and its host galaxy metallicity 
is not surprising in the context of the collapsar model, which 
requires rapidly-rotating stars that lack a H/He envelope (in 
order to be in accord with supernova observations). These 
requirements can be satisfied by a metal-poor progenitor star. 
Introducing a metallicity bias leads to an accelerated evolution 
of the cosmic GRB rate density (relative to the SFH), which 
allows for a better fit to recent Swift data. In models that at- 
tribute ultrahigh energy cosmic rays to gamma-ray bursts, this 
evolution provides the history of UHECR production in the 
universe. If GRBs are to account for the observed cosmic -ray 
spectrum, they must generate a flux of cosmogenic neutrinos. 

A broad approach which utilizes all available means of 
observation will continue to be of great utility in unveiling 
the mysteries of gamma-ray bursts, particularly a relation to 
UHECR. Further observations of GRB host galaxies, along 
with those of core-collapse supernovae (particularly Type 
Ib/Ic), will provide invaluable information concerning GRB 
progenitors. A systematic study of host metallicities would 
be important for firmly establishing the GRB -metallicity anti- 
correlation, which, along with more accurate measurements of 
cosmic metallicity evolution 115011 . would allow for improved 
GRB rate-related calculations. 

Direct observations of neutrinos from GRBs (e.g., by Ice- 
Cube [40] or a km 3 Mediterranean detector B41I1 ) would con- 
firm the acceleration of protons to high (though not necessar- 
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ily ultrahigh) energies. In order to produce a detectable signal, 
a particularly strong burst, with measured gamma-ray fluence 
> 3 x 10: 4 erg cm" 2 (e.g., the recent GRBs 060928 JH and 
061007 B85I1 . which were observed to have fluences very near 
this threshold) is required j62ll . While the needed neutrino 
telescopes (which may also reveal the sources of Galactic cos- 
mic rays Km) are still under construction, future observations 
by Swift 115611 should reveal similar bursts, in addition to deter- 
mining the redshift-dependent GRB rate with unprecedented 
precision. 

Cosmogenic neutrinos present a unique tool to examine the 
GRB-as-proton accelerator conjecture. While the cosmic-ray 
spectrum that will be measured by Auger will allow for fur- 
ther assessments of the viability of prospective source models, 
the combined measurements of neutrinos and UHECR would 
break degeneracies between the various models, as discussed 
by Seckel and Stanev |3j]- It is well-known that strong source 
evolution can lead to observable cosmogenic neutrino signals; 
however, a physically motivated model was lacking. The en- 
hanced cosmogenic neutrino fluxes expected to result from 
strong GRB evolution will allow for this model to be tested 
in a novel fashion. We have produced this strong evolution 
naturally through metallicity dependence, however, our gen- 
eral result is also applicable to any GRB model that attempts 
to explain the observed Swift redshift distribution with addi- 
tional evolution. 

The sensitivity afforded by ANITA will allow for near-term 
examination which may either affirm or, if the expected flux 
is not found, place substantial constraints upon the model 
parameters. The lower detection threshold achieved by an 
ARIANNA-like detector will allow for a realistic opportu- 
nity to discriminate between evolution models. Measurements 
of fluxes consistent with that expected from GRB evolution 
would provide compelling evidence for this model. As the 
relatively young field of particle astrophysics continues to 
progress, it is quite an exciting possibility that metals, once 
relegated to the realm of pure astronomy, may, in fact, hold 
the key to the production of the highest energy particles in the 
universe. 
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APPENDIX A: GRB REDSHIFT DISTRD3UTION 

In calculating the redshift distribution of gamma-ray bursts 
expected to be observable by Swift, we follow the model 
of Le and Dermer Olll . using similar, but somewhat simpli- 
fied, notation. This requires the assumption of a character- 
istic GRB gamma-ray energy output (£y) and a duration in 
the GRB rest frame (8t), additionally assuming flat temporal 
and spectral burst profiles. Each burst is assigned a jet open- 
ing angle, 9, which is selected from a power law distribution 
of the form go(l — cos 9) s , between an assumed cos max and 



cos0 m i n , normalized with go = (1 — cos0 max ) 

(l-COSftnU,) 1 "*]. 

The observed energy flux from a GRB at a luminosity dis 
tRnce, d((z), then takes the form 



1+1 _ 



tance 



By /St 



\%d)(\ -cos0)A/, 



(Al) 



where Xb — 5 is a bolometric correction factor, which accounts 
for the fraction of the burst spectrum within the energy band of 
the GRB detector. Solving this equation for the detector trig- 
ger sensitivity threshold, /d e t, yields the beaming angle, 0det. 
required for detection. In effect, by assuming a constant £y we 
will use the opening angle as a proxy for the isotropic equiva- 
lent energy (Ei w ) often used in GRB studies (e.g., Ref. [87]). 

For a comoving GRB rate, n c (z), the detectable GRB rate is 
reduced to n c (z)(l — cos 6) by beaming through angle, 9. We 
then find the detectable GRB fraction as 



©(z) = / (1 -cose)go(l -cos0)' s ,icos0 

JcosMin[9 max ,e det ] 



1 +s (1 -cosMin[0 max ,0 det ]) 2+4 '-(l -cos0 min ) 



\2+s 



2 + s 



(1 -cos0 max ) 1+1 - (1 -cos0 min ) 



L+s 



(A2) 

We then compute the number of detected bursts per unit red- 
shift as 



dN l CM\ n °® dV c. 

&{z) - ■ _ —f-dz, 



dzdQ. 4n 



l+z dz 



(A3) 



where 1/(1 + z) takes into account cosmic time dilation, and 
the comoving volume element is given as 



dV c _ \%d\c 



dz {\+z)dz/dt 



(A4) 



Substitution of Eqs. dA2b and ( lA4b into Eq. ( IA3b agrees 
with the formulation presented in Ref. f3~Ill . 



APPENDIX B: COSMOGENIC NEUTRINO SPECTRUM 

We describe a simple, analytic method to compute the cos- 
mogenic neutrino spectrum observed at Earth, based upon the 
assumption that each daughter neutrino resulting from a pho- 
topion interaction receives ~ 1 /20 of the parent proton energy. 
That is, a neutrino detected with energy, E v , was produced 
with energy, E' v = (1 + z)E v , at the source, from a parent pro- 
ton which had energy, E' p = 2QE' V — 20 (1 +z)E v . We predict 
the observed (z = 0) neutrino spectrum, dN v /dE v , by integrat- 
ing over the contributions to the flux from redshifts up to Zmax 
as 



dN v 
dE v 



c 

An 



dN v dE' 



dE' v dE v dz/dt 



dz 



(Bl) 



where dE' v j dE v = (l+z) accounts for the fact that neutrinos 
are observed in a narrower energy bin then they were origi- 
nally produced (due to redshifting). Evolution in the density 
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of sources is accounted for through W{z). The neutrino spec- dN v dE' v _ (pi\<y\(]j_ \ p'l cri\ 

trum at production, dN v /dE' v , can be directly related to the l¥ v ~dE~ v 2X3 ^pl-^ \\ l +z > x E p\ 

cosmic-ray injection spectrum <p{E' p ) (per unit comoving vol- 
ume per unit energy, per unit time), through the conservation 
of the transferred energy, as 

,iy ^jji with as defined in Section HVI Substitution of Eq. ( |B3t into 

E ' V W = ^ (P( ^ ) ^[ (1+Z)X ^])^ (B2) Eq. (EB yields Eq. ©. 
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